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Motivated by the fact that cahbrated light curves of Type la supernovae (SNe la) have become 
a major tool to determine the expansion history of the Universe, considerable attention has been 
given to, both, observations and models of these events over the past 15 years. Here, we summarize 
new observational constraints, address recent progress in modeling Type la supernovae by means 
of three-dimensional hydrodynamic simulations, and discuss several of the still open questions. It 
will be be shown that the new models have considerable predictive power which allows us to study 
observable properties such as light curves and spectra without adjustable non-physical parameters. 
This is a necessary requisite to improve our understanding of the explosion mechanism and to settle 
the question of the applicability of SNe la as distance indicators for cosmology. We explore the 
capabilities of the models by comparing them with observations and we show how such models can 
be applied to study the origin of the diversity of SNe la. 
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I. INTRODUCTION 

Today, Type la supernovae (SNe la for short) play a 
somewhat ambiguous role in astrophysics. On the one 
hand, it is their relative homogeneity that caused their 
use as distance indicators in observational cosmology. On 
the other hand, this evoked an enormous interest result- 
ing in a rather comprehensive observational survey of 
SNe la that over the last decade clearly revealed sub- 
classes with diverging properties and variability among 
these objects. While the notion of homogeneity inspired 
the model of SNe la being explosions of Chandrasekhar- 
mass white dwarfs (WDs), the newly discovered hetero- 
geneity of the class suggests multiple progenitors and/or 
explosion mechanisms. 

After the pioneering work by Arnett [I] numerical sim- 
ulations have been instrumental in modeling supernovae. 
Until the 1990s this approach was restricted to one spatial 
dimension which prevented a realistic treatment of the 
multi-dimensional burning mechanism in these objects. 
However, parametrized models of that time, notably the 
W7 model of Nomoto et al. [5], still set a standard in the 
field and are widely used in the interpretation of obser- 
vational data. In the 1990s the first multi-dimensional 
SN la simulations emerged. Together with earlier work 
on one-dimensional models, they are reviewed by Hille- 
brandt and Niemeyer [3]. Here, we report on develop- 
ments in the last decade, focusing on work associated 
with the supernova group at the Max-Planck-Institute 
for Astrophysics, Garching, but putting it into context 
with other work. 

While the modeling of the explosion physics has made 
substantial progress (in particular with the introduction 
of multi-dimensional simulations), the question of the 
progenitor system of SNe la remains a fundamental prob- 
lem. There is wide consensus that these events are due 
to thermonuclear explosions of WDs P| , most likely com- 



posed of carbon and oxygen. This was recently confirmed 
by Nugent et al. [S] and Bloom et al. [S] who on the basis 
of early time observations concluded that the exploding 
object in SN la 201 Ife was a compact star. The question 
of how it reaches an explosive state, however, is more 
complicated. As single WDs are unconditionally stable, 
some kind of interaction with another star is necessary to 
explain the supernova. Unfortunately, attempts to iden- 
tify this star beyond doubt have failed so far. Collisions 
with compact objects in globular clusters [7]-(5] lead to 
atypical events or fall short of explaining the SN la rate. 
Thus, although such events may occur in Nature, the 
bulk of SNe la is more likely to be associated with stellar 
binaries. The nature of the binary companion, however, 
is still unclear. Traditionally, two classes of potential bi- 
nary progenitor systems have been distinguished - the 
single-degenerate progenitor channel, in which the com- 
panion is a normal star, and the double-degenerate chan- 
nel with two WDs interacting and merging. At present 
it is unclear whether one of these possibilities is exclu- 
sively realized in Nature or whether both contribute to 
the class of SNe la. 

This problem has been approached from different per- 
spectives. Observational data becomes increasingly con- 
straining for the physical mechanism of SNe la and a brief 
overview of the current status is given in Sect. [H] In ad- 
dition, the rate at which SNe la occur and the distribu- 
tion of delay times between formation of the progenitor 
systems and supernova explosions can help to identify 
the dominant progenitor channel(s). These data can be 
compared with predictions from binary population syn- 
thesis calculations. We discuss this approach and recent 
results in detail in Sect. |III[ Another possibility is to fol- 
low different explosion scenarios in hydrodynamic sim- 
ulations. Combined with radiative transfer calculations 
these predict observables that can be directly compared 
to SN la observations. Over the last decade substan- 
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tial progress in this approach was possible due to fully 
multi-dimensional treatment that allows to reduce the 
free parameters involved in describing the physics and 
thus improve the predictivity of the models. We dis- 
cuss recent results in Sect. |IV[ where the modeling ap- 
proaches and the implementation in numerical simula- 
tions are briefly outlined followed by the presentation of 
models that potentially can account for normal SNe la 



(Sect. IV B I, while other models seem to reproduce pecu- 



liar subclasses (Sect. IV C I. We emphasize that this is a 



way of presenting our models and discussing the results. 
However, the assignment of models to different subclasses 
is not necessarily unique, but it is chosen here to point 
out the possibility to model a wide variety of SNe la when 
considering different progenitor scenarios and explosion 
mechanisms. 



II. OBSERVATIONS 

The efforts to systematically obtain observational data 
of SNe la have gained tremendous momentum during the 
past 15 years. This is primarily a result of their un- 
equaled potential to act as 'standardizable' candles for 
the measurement of the cosmological expansion rate and 
its variation with look-back time PHH^ (see also Goobar 
and Leibundgut [22] for a recent review). The discovery 
that the Universe entered into a phase of accelerated ex- 
pansion at a redshift of around 0.5, due to the action of 
some unknown form of 'dark energy', was awarded with 
the Nobel Prize in Physics in 2011 to Saul Perlmutter, 
Adam Riess, and Brian Schmidt. 

For theorists, this development presents both a chal- 
lenge, to help to understand the correlations among the 
observables, and an opportunity, to use the wealth of new 
data to constrain the zoo of existing explosion models. 
There exist a number of excellent reviews about SNe la 
observations in general [23H25] . their spectral properties 
[26] [27] . and photometry in the IR and optical bands 
[^5H5D] . Here, we highlight those aspects of SN la ob- 
servations that most directly influence theoretical model 
building at the current time. 



A. General properties 

The classification of SNe la is based on spectroscopic 
features: the absence of hydrogen absorption lines, dis- 
tinguishing them from Type II supernovae, and the pres- 
ence of strong silicon lines in their early and maximum- 
light spectra, classifying them as Type la's [ST] . 

The spectral properties, absolute magnitudes, and 
light-curve shapes of the majority of SNe la are remark- 
ably similar, exhibiting only small spectroscopic and pho- 
tometric differences [12]. It was believed until recently 
that approximately 85% of all observed events belong to 
this class of 'normal' 32J SNe la, represented for example 
by SN 1972E, SN 1994D, or SN 2005cf. However, recent 



studies show that the peculiarity rate can be as high as 
30% as suggested for instance by Li et al. [33] . 

The optical spectra of normal SNe la contain neutral 
and singly-ionized lines of Si, Ca, Mg, S, and O at max- 
imum light, indicating that the outer layers of the ejecta 
are mainly composed of intermediate mass elements [26] . 
Permitted Fe ii lines dominate the spectra roughly two 
weeks after maximum when the photosphere begins to 
penetrate Fe-rich ejecta [371 El]- In the nebular phase, 
beginning approximately one month after peak bright- 
ness, forbidden Fell, Feiii, and Com emission lines be- 
come the dominant spectral features. Some Can remains 
observable in absorption even at late times [5^ . The de- 
crease of Co lines and the relative intensity of Co ill and 
Feiii give evidence that the light curve tail is powered 
by radioactive decay of ^^Co ^(see also Truran et al. 
[26], Colgate and McKee [37]). 

The early spectra can be explained by resonant scatter- 
ing of a thermal continuum with P Cygni-profiles whose 
absorption component is blue-shifted according to ejecta 
velocities of up to about 25,000 km s~^, rapidly decreas- 
ing with time. Different lines have different expansion 
velocities 38^30], suggesting a layered structure of the 
explosion products. 

Photometrically, SNe la rise to maximum light in a 
period of approximately 18 to 20 days [TOffi] reaching 
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^19.30 ±0.03 + 5 log(iJo/60) (1) 



with a dispersion of cta/ < 0.3 [45]. It is followed by a 
first rapid decline of about three magnitudes in a matter 
of one month. Later, the light curve tail falls off in an 
exponential manner at a rate of approximately one mag- 
nitude per month. In the I- and near-infrared bands, 
normal SNe la rise to a second maximum approximately 
20 days after the first one [H]. Typical ^^Ni masses in- 
ferred from their bolometric light curves are in the range 
from 0.3 to 0.9 for normal SNe la (e.g., Stritzinger 
et al. [46]). 

It is especially interesting that the two most abundant 
elements in the universe, hydrogen and helium, so far 
have not been unambiguously detected in the spectra of 
normal SNe la [Ml HZ] (but see Meikle et al. [35] and 
Mazzali and Lucy [H] for a possible identification of He, 
and Hamuy et al. [4^ and Dilday et al. [50] for an iden- 
tification of H in individual interacting peculiar objects). 
Also, there are no indications yet for radio emission [51] . 
including the rather nearby supernova SN 201 Ife [531155] . 



B. Diversity and correlations 

Early suggestions [M] [55] that the existing inhomo- 
geneities among SN la observables are strongly inter- 
correlated are now established beyond doubt [321 US]- 
Branch [12] summarizes the correlations between spec- 
troscopic line strengths, ejecta velocities, colors, peak 
absolute magnitudes, and light curve shapes that were 
known at that time. Roughly speaking, SNe la appear 
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FIG. 1. Observational diversity of SNe la in B-band decline 
rate Ami5(_B) and B-band peak absolute magnitude. Normal 
SNe la (shown in grey, data taken from Hicken et al. |57]) fol- 
low the Phillips relation !58]. 1991bg-like SNe (shown in red, 
data taken from Taubenberger et al. j^), and 2002cx-like 
SNe (shown in green, data taken from Phillips et al. [60]) are 
subluminous with respect to the Phillips relation. Superlumi- 
nous SNe la (shown in cyan, data from Taubenberger et al. 
|61| ) are almost one magnitude brighter in B-band than nor- 
mal SNe la with a comparable B-band decline rate. 



to be arrangeable in a one-parameter sequence according 
to explosion strength, wherein the weaker explosions are 
less luminous, redder, and have a faster declining light 
curve and slower ejecta velocities than the more ener- 
getic events. Based on these findings Mazzali et al. [56] 
argue that a single explosion scenario, possibly a delayed 
detonation (see Sect. IVBl), may explain most SNe la. 
However, more recent (and better) data challenge this 
conclusion, as will be discussed below. 

The relation between the width of the light curve 
around maximum and the peak brightness (brighter su- 
pernovae decline more slowly) is the most prominent of 
all correlations (Fig.[l] Phillips [TO, Pskovskii [53]). Pa- 
rameterized either by the decline rate Amis [lOlllS], a 
'stretch parameter' |62| . or a multi-parameter nonlinear 
fit in multiple colors jllj . it was used to renormalize the 
peak magnitudes of a variety of observed events, sub- 
stantially reducing the dispersion of absolute brightnesses 
(see, e.g., Leibundgut and Goobar and Leibundgut 
[22j for recent reviews). This correction procedure is a 
central ingredient of all current cosmological surveys that 
use SNe la as distance indicators [63H68] . 

However, there are supernovae, classified as Type la, 
which violate this correlation. SN 1991bg, SN 1992K, 
SN 1999by and SN 2005bl are well-studied examples 



for red, fast, and subluminous supernovae with a typi- 
cal Ami5 value of about 1.8 and B-band peak absolute 
magnitudes around -17, roughly one magnitude fainter 
than their 'normal' counterparts [551 [S51 -f75] . Their V , /, 
and i?-band light curves decline unusually quickly, skip- 
ping the second maximum in /, and their spectra show a 
high abundance of intermediate mass elements (including 
Ti 11) with low expansion velocities but only little iron. 
Models for the nebular spectra and light curve of SN 
1991bg consistently imply that the total mass of ^^Ni in 
the ejecta was very low (~ 0.07Mq[73]), a typical value 
for this class being ^ 0.1 M©. In addition, there is also 
evidence for unburnt C and O in their early spectra, in 
contrast to normal SNe la. These 'subluminous' explo- 
sions make up for about 15% (or more) of all SNe la 

The prototype of a second group of subluminous SNe la 
is SN 2002CX [75l |76]. Here, again, the mass of ^^Ni, 
as estimated from 'Arnett's rule' [77], is low, around 
0.2 M0 only. The spectra show narrow lines, indicating 
low ejecta velocity and low kinetic energy. Other super- 
novae belonging to this class include SN 2005hk ^60^ .78] , 
SN2008ge [70], and SN 2009ku [SO]- According to Li 
et al. [33] they contribute about 5% of all SNe la. Even 
300 days after the explosion, the ejecta of members of 
this group are not transparent, but show emission from 
a narrow region in velocity space (less than 1000 km s~^ 

m)- 

Finally, transients even fainter than 1991bg-like SNe 
have been observed, SN 2005E [51] or SN 2005cz [M] 
being examples. They are Ca-rich fast decliners, their 
spectra resemble more SNe lb than SNe la, i.e., they 
show He but little O and Si in their early-time spectra, 
and their decline rates are similar to those of SNe Ic. 
They are found in old stellar populations, however, and 
the discussion is open whether they are thermonuclear 
explosions or core-collapse supernovae [83H85] . 

At the other end of the luminosity function, SN 1991T 
is often mentioned as a striking representative of bright, 
energetic events with broad light curves [69l [86U89j . 
Rather than the expected Sin and Can, its early spec- 
trum displayed high-excitation lines of Fe ni but returned 
to normal a few months after maximum light. But re- 
cently other SNe la were found which are even more lu- 
minous than SN 1991T, with decline rates that put them 
well above the Phillips relation by almost one magnitude 
in the i3-band, prototypical examples being SN 2006gz 
and SN 2009dc [SH ^MM- By now, seven objects that 
may belong to this subclass have been discovered and 
they may contribute up to about 9% of all SNe la j33j . 
In addition to their high luminosity, 2 to 3 times higher 
than normal SNe la, they are characterized by a slow 
decline {/!\mi^{B) ~ 0.8), a long rise time (> 23 days), 
low ejecta velocities, and prominent C n absorption fea- 
tures, while other properties of their early-time spectra 
are similar to those of normal SNe la. If the luminosity at 
peak would come exclusively from the decay of ^^Ni the 
Ni-mass of SN 2009dc would be around 1.5 to 1.8 M© 
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[5T1 1^ . exceeding the Chandrasekhar mass (see, how- 
ever, Hachinger et al. |96j for an alternative scenario). 
The various sub-classes are illustrated in Fig. [T] 

From early on, peculiar events like SN 1991T and SN 
1991bg were suggested to belong to different subgroups 
of SNe la than the normal majority, created by different 
explosion mechanisms [BSl [ZH [57] although the overall 
SN la luminosity function seems to be rather smooth, 
with a shallow increase from an absolute i?-band magni- 
tude of -17 to -19, followed by a steep decline to -19.5 [33] 
(thus leaving out 09dc-like events), indicating that 'nor- 
mal' SNe la are essentially the brightest, with Ni masses 
around 0.6 while the full class may contain a large 
number of undetected subluminous SNe la. 



C. SNe la and their host galaxies 

There is mounting evidence that SN la observables 
are correlated with their host stellar population [T^ and 
there are recent investigations demonstrating the depen- 
dence of supernova properties on global characteristics 
of their hosts [55Vll02| . For instance, SNe la in red or 
early-type galaxies show, on average, slower ejecta veloc- 
ities, faster light curves, and are dimmer by ~ 0.2 to 0.3 
mag than those in blue or late-type star-forming galax- 
ies [ISlESllMllinilinil- Moreover, SNe la seem to have 
lower ejecta velocity in high-mass host galaxies |105j . On 
the other hand side, SNe la at low and high redshift seem 
to have similar spectral evolution |1061 1107j . 

The SN la rate per unit stellar mass is nearly a fac- 
tor of 20 higher in late-type galaxies than in early-type 
ones and depends inversely on the host galaxy's mass 
[108H110) . The rate seems to be lower in galaxy bulges 
than in spiral arms. These findings indicate that there 
might be a population of progenitors with large delay 
time |1111 I112j . Also, the outer regions of spirals ap- 
pear to give rise to similarly dim SNe la as ellipticals 
whereas the inner regions harbor a wider variety of ex- 
plosion strengths |113j . 



D. Summary: observational requirements for 
explosion models 

To summarize the main observational constraints, any 
viable scenario for the SN la explosion mechanism has 
to satisfy the following (necessary but probably not suf- 
ficient) requirements; 

1. Agreement of the ejecta composition and velocity 
with observed spectra and light curves. In gen- 
eral, the explosion must be sufficiently powerful 
(i.e., produce enough ^^Ni) and produce a substan- 
tial amount of high- velocity intermediate mass ele- 
ments in the outer layers. Furthermore, the isotopic 
abundances of "normal" SNe la must not deviate 
significantly from those found in the solar system. 



2. Robustness of the explosion mechanism. In order 
to account for the homogeneity of normal SNe la, 
the standard model should not give rise to widely 
different outcomes depending on the fine-tuning of 
model parameters or initial conditions. 

3. Intrinsic variability. While the basic model should 
be robust with respect to small fluctuations, it must 
contain at least one parameter that can plausi- 
bly account for the observed sequence of explosion 
strengths. However, this could in principle also be 
achieved by allowing for different progenitor chan- 
nels. 

4. Correlation with progenitor system. The explosion 
strength parameter must be causally connected 
with the state of the progenitor WD in order to 
explain the observed variations as a function of the 
host stellar population. Moreover, there must be a 
sufficient number of progenitor systems such that 
the rate and delay-time constraints are matched. 

III. PROGENITORS 

While it is widely accepted that SNe la originate from 
explosions of WDs that approach critical conditions such 
that burning can proceed explosively (cf. Sect. |l]), the 
manner in which these conditions are achieved remains 
uncertain. Almost certainly the WD gains matter from a 
nearby stellar companion. Until very recently, the stan- 
dard paradigm was the following: SNe la originate from 
probably one, possibly two different formation channels 
which enable the WD to reach critical conditions neces- 
sary for a thermonuclear explosion to occur (see following 
subsections). 

However, recent observations of SNe la have brought to 
light the (previously shrouded) highly diverse nature of 
these objects (see Sect.|ll|. When one considers all of the 
necessary criteria that a progenitor model must satisfy in 
order to be seen as a viable progenitor candidate (robust 
explosion mechanism, ejecta stratification, velocities and 
nucleosynthesis, characteristic peak luminosity and light- 
curve shape, absolute birth rates and delay times), it is 
evident that reconciliation of the entire range of observed 
characteristics of SNe la with a single progenitor scenario 
is improbable. If more than one progenitor scenarios are 
contributing to the observed population of SNe la - which 
is currently the favoured view |114j - it is still unclear as 
to which progenitor scenario(s) dominate(s). 

Binary population synthesis models have been used for 
a few decades now to estimate relative (and absolute) 
birthrates of various binary formation channels that can 
lead to SNe la (see |115^ for one of the most well-known 
early studies). A powerful feature of population synthe- 
sis models is that one is able to easily compute the delay 
times of SNe la, which puts strict limits on the system 
age, thus ruling out certain theoretical progenitor sce- 
narios. In addition, the models enable one to reconstruct 
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the entire evolutionary history for all binaries of interest, 
which is critical for uncovering evolutionary phases (e.g. 
mass transfer episodes) that might give rise to observa- 
tional features which could kill or confirm a given model 
(see e.g. Ruiter et al. ^llGj in which the StarTrack [117] 
binary evolution synthesis code is used). 

In terms of SN la rates and delay times, the results 
from different population synthesis codes are found to 
vary quite a lot in some cases and agree fairly well in oth- 
ers (Nelemans et al. |118j and references therein). This 
is primarily due to the rather uncertain nature of mass 
transfer and accretion in close binary stars which leads 
to differing assumptions for the input physics in the var- 
ious codes. In particular, for progenitors which undergo 
(quasi)stable mass exchange, one must decide how the 
(donor) mass transfer rate is approximated in the given 
binary and correspondingly, how the (accretor) retention 
efficiencies are computed (see, e.g., |119l ll2d] V Further 
still, the manner in which matter is lost from the binary 
(carrying away with it angular momentum) will also have 
an effect on the orbital behaviour and subsequent binary 
evolution (121| . 

The progenitor problem is still unsolved, though as 
previously mentioned, it seems likely that at least two 
progenitor scenarios (and possibly more explosion mech- 
anisms) are required in order to explain the observed SN 
la rate and delay time distribution [122| . In the follow- 
ing sub-sections we review the most promising progenitor 
scenarios (e.g. formation channels) which are thought to 
lead to SNe la. 



A. Single-degenerate (Chandrasekhar-mass white 
dwarf) scenario 

1. Hydrogen-burning donors 

Often called the single degenerate (SD) scenario |123j . 
SD systems are detected in Nature (e.g. RS Oph, [TM]) 
and were, in the past, widely thought to be the most 
promising SN la progenitors. In this scenario, the com- 
panion star is a main sequence or giant-like star (possibly 
a helium-burning star; see Sect. IIIA2I that is overfill- 



ing its Roche-lobe, transferring matter through the in- 
ner Lagrange point in a stable manner to the companion 
carbon-oxygen (CO) WD. If the mass transfer proceeds 
within a certain range of rates (for example |125p . the 
donor material is accreted in a stable fashion leading to 
efficient hydrogen-burning (mass accumulation) on the 
WD, thereby increasing its (central) density. When the 
density in the center of a CO WD becomes high enough 
the carbon in the WD starts to burn (see Sect. IV) which 



eventually leads to a thermonuclear explosion, obliterat- 
ing the WD and possibly imparting a significant kick 
on the companion star (Sect. |l]). This critical density 
when carbon-burning can start is usually attained when 
the WD approaches a critical mass - the Chandrasekhar 
limit. 



A typical formation pathway leading to the SD sce- 
nario involves an episode of unstable mass transfer fol- 
lowed by an episode of stable mass transfer at a later 
stage. A more specific example (e.g. jl26l ) is as follows: 
the initially more massive star (the primary) first loses its 
hydrogen-rich envelope on the asymptotic giant branch 
(AGB) when it fills its Roche-lobe and mass transfer is 
dynamically unstable. This results in a common enve- 
lope, which serves to bring the two stars to a smaller or- 
bital separation |127j . The post-common envelope binary 
comprises a (newly formed) CO WD and a (likely still on 
the main sequence) companion. At some later stage, the 
companion then fills its Roche-lobe (either while on the 
main sequence or as an evolved star), only this time mass 
transfer is stable, and the CO WD grows in mass until it 
approaches the Chandrasekhar limit. 



2. Helium-burning donors 

It is also possible that the WD may reach the Chan- 
drasekhar mass, by accreting from a helium-burning star 
donor rather than a hydrogen-rich donor (e.g. |128| ). 
Such formation channels are expected to be rare, and 
these progenitors have shorter evolutionary timescales 
than the 'canonical' SD scenario due to the larger 
zero-age main sequence (ZAMS) mass of the secondary 
[1261 11291 1130j . Since they might also harbour distinctly 
different physical (observable) properties, they have been 
considered their own class by some authors (e.g., helium- 
rich (HeR) scenario, |131) ). 



B. Double-detonation (sub-Chandrasekhar-mass 
white dwarf) scenario 

Another progenitor channel which has recently 
re-gained popularity among the community is the 
sub-Chandrasekhar-mass scenario in which a sub- 
Chandrasekhar-mass CO WD accretes stably from a 
companion and never reaches the Chandrasekhar limit 
before exploding [129] . Depending on the assumed mass 
transfer/accretion rates and the mass of the CO WD, 
the WD is thought to be able to accumulate (rather 
than burn) a layer of helium which may detonate un- 
der the right physical conditions [1321 1133j . This shell- 
detonation - if realized - likely triggers a second deto- 
nation in the sub-Chandrasekhar-mass WD, leading to a 
SN la ('double-detonation' scenario, see Sect. IVB2). 

Such a scenario was investigated from a population 
synthesis standpoint by Tutukov and Yungelson |134j in 
context of helium-rich donors, and by Yungelson et al. 
jl35| in context of of symbiotic systems. However, forma- 
tion channels leading to the double-detonation scenario 
via accretion from a hydrogen-rich companion are chal- 
lenged by the ability of the WD to efficiently accrete (and 
stably burn) hydrogen and helium |136[ 1137) . Nonethe- 
less, a double-detonation progenitor scenario might be 
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readily realized in Nature from sub-Chandrasekhar-mass 
WDs accreting from helium-burning stars or helium-rich 
degenerate (or semi-degenerate) stars [1341 1138] . 

In considering all potential helium-rich donors trans- 
ferring mass to sub-Chandrasekhar-mass CO WDs, 
Ruiter et al. [131j investigated the double-detonation sce- 
nario. In that work, two characteristic evolutionary chan- 
nels were found (note that the authors only considered 
SNe la to arise from systems where the total WD mass 
was > 0.9 Mq). The typical formation channels were as 
follows: 

Helium-burning star donors. This formation channel in- 
volves two cases of unstable mass transfer (common en- 
velopes) followed by a later stage of stable mass transfer 
from the secondary to the CO WD. The initial primary 
star fills its Roche-lobe (unstably) while on the AGB 
when the companion is still on the main sequence, re- 
sulting in a CO WD-main sequence binary. A second 
common envelope ensues when the secondary star - now 
evolved - fills its Roche-lobe. The post-common envelope 
binary consists of a CO WD (left over from the first com- 
mon envelope) and a stripped core of a giant; a compact 
naked helium-burning star. Since the stars are already 
on a fairly close orbit (due to two common envelopes), 
it does not take long for the stars to be brought into 
contact. The naked helium-burning star then fills its 
Roche-lobe and transfers matter stably to the primary 
WD, until the onset of the double-detonation. 
(Semi-)degenerate star donors. A number of evolution- 
ary pathways can lead to the formation of such a progeni- 
tor, but the most common path also involves two common 
envelopes (primary on the AGB, then secondary on the 
giant branch), followed by a phase of stable mass trans- 
fer. However, in this case, the secondary's ZAMS mass 
is smaller than that of helium-burning star donor case. 
Thus, when the secondary loses its H-rich envelope in 
the second common envelope, a degenerate (non-burning) 
naked helium core of a giant is left behind. Once contact 
is achieved, the 'helium-rich WD' transfers matter stably 
to the CO WD until the onset of the double-detonation. 



C. Double degenerate mergers 

Another scenario which might readily lead to SNe la 
is the merger of two CO WDs where the total mass ex- 
ceeds the Chandrasekhar limit ('double degenerate' (DD) 
scenario, |139j ). Along with the SD scenario, the DD 
scenario has been a leading progenitor candidate model. 
The reason is owed partially (but not only) to the theoret- 
ical birth rate calculations - for which it historically does 
the best of any progenitor scenario. A number of pop- 
ulation synthesis studies over the last few decades (e.g. 
[115) ) have demonstrated that mergers of CO WDs with 
a total mass exceeding the Chandrasekhar limit might 
be frequent enough to account for Galactic SN la rates, 
depending on e.g. the adopted prescriptions for common 
envelope evolution |126U140] . Still, when considering cos- 



mological SN la rates as a function of delay time, the DD 
model scenario, like other scenarios, often falls short of 
the observationally-recovered rates by at least a factor 
of a few |131l 11411 I142j . However, some of the most 
recently-measured delay time distributions indicate that 
DD merger rates might indeed be frequent enough to 
account for the bulk of SNe la at least in some stellar 



populations ([Iinillll], see Sect. jlllE l. 

There are a number of progenitor pathways that can 
lead to a CO-CO WD binary. A typical one would be 
the following |126) : The primary star fills its Roche-lobe 
when it is slightly evolved, and mass transfer is stable 
to the companion (a second stable phase of Roche-lobe 
overflow may follow a bit later, when the primary is an 
evolved helium star) . The primary star then evolves into 
a CO WD. The secondary star fills its Roche-lobe when 
it is an evolved star but mass transfer is unstable to 
the WD, and a common envelope ensues. The post- 
common envelope binary consists of a CO WD and a 
naked helium-burning star. Following this, a final phase 
of stable mass transfer may occur whereby the slightly 
evolved helium star (secondary) transfers matter to the 
primary WD. Such a final phase of mass transfer was 
found to be important in explaining the peak-brightness 
distribution of SNe la (see Sect. jlVB 4] ). Once the WDs 
reach contact and the larger (less massive) one fills its 
Roche-lobe, mass transfer must be unstable for a merger 
to occur. This is the likely outcome for double CO WDs 
given their typical mass ratios (see [143j : see also Too- 
nen et al. (141) for specific examples of DD formation 
channels) . 

Despite favourable theoretical rate predictions, the DD 
scenario has received a lot of criticism over the years. 
Earlier calculations predicted that the likelihood of ther- 
monuclear explosion in a double CO WD merger is rather 
unlikely. Such a merger (Mtot ^ Mch) was thought to 
lead to disruption of the secondary WD which is then ac- 
creted onto the primary. The accretion would not lead to 
central burning but rather burning in the outer layers of 
the WD, where densities are lower, and the accreting CO 
WD would transform into an oxygen-neon-magnesium 
WD |144| . For WDs of such composition, electron cap- 
tures become important at high central densities, and as 
the WD approaches the Chandrasekhar mass it collapses 
to form a neutron star in an accretion-induced collapse 
(AIC, mSj). Even if DD mergers do lead to SNe la, it 
has been reiterated by some groups that, assuming the 
estimated birthrates from population synthesis calcula- 
tions are correct within a factor of a few, all CO WD 
mergers (even those with Mtot < 1-4 Mq) must lead to 
SNe la in order to match the observed rates |1461 1147j . 
It is unclear whether lower-mass mergers lead to ther- 
monuclear explosions, let alone if such explosions would 
produce enough ^^Ni. 

Depending on the configuration of the binary system - 
in particular the mass ratio - the merger may be some- 
what quiescent as described above, or it may be violent 
enough such that a prompt detonation in the primary 
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WD will occur. These 'violent mergers' are robustly 
found to lead to a thermonuclear explosion, and they 
are described in Sect. IIV B 31 



D. Other possible scenarios 

Other possible formation channels leading to SNe la 
have been postulated in the literature. For example: 
a potential scenario involves the merger of a CO WD 
and the core of an AGB star during a common envelope 
event |148l 1149] . Such events are expected to readily oc- 
cur, however it is unclear whether such 'core-degenerate' 
mergers would lead to an immediate (or delaved. jlSOj ) 
explosion that exhibits observational signatures which 
match those of SNe la. 



E. Constraining progenitor models: delay times 
and rates 



The delay time distribution (DTD) is the distribution 
of times in which SNe la explode following a (hypotheti- 
cal) burst of star formation. Knowing the DTD gives the 
age of the progenitor, which places strong constraints on 
the different proposed progenitor scenarios. If the SN la 
rate is known in addition, then it becomes possible to 
rule out theoretical formation channels. 

Calculation of observationally-recovered DTDs in- 
volves many assumptions, the most important being the 
assumed star formation history of the supernova's host 
galaxy or local stellar population, for which several tech- 
niques have been employed [1101 I151j . There are two 
emerging facts in the literature: i) there is a population 
of 'prompt' SNe la which have delay times ^500 Myr and 
seem to comprise a significant fraction of all SNe la, and 
a) there are SNe la which are 'delayed', the seemingly- 
continuous DTD spanning up to a Hubble time, with a 
characteristic (cosmic) DTD beyond > 400 Myr that fol- 
lows a power-law shape |142| . Such a power-law 
(~t~^) is expected if the dominating timescale leading 
to SNe la is set by gravitational radiation, as is the case 
for DD mergers jl52j . 

With binary population synthesis models, the entire 
evolutionary history of each binary is followed, so the 
DTD is easily determined for all potential SN la pro- 
genitors. In Figure [2] we show mass-normalized the- 
oretical DTDs for the following progenitor scenarios: 
SD (hydrogen-burning and helium-burning donors); a 
sub-class of the DD, whereby the mass of the primary 
WD must be > 0.8 M0and there are additional restric- 
tions on the mass ratio, (cf. Sect. IVB3I; and double- 
detonation progenitors involving both helium-burning 
star and helium-rich WD donors [1161 11261 1131j . Along- 
side our theoretical DTDs we show the most recent DTDs 
from (red squares, Maoz et al. |153j l and the best-fit 
DTD from observations as described in Graur and Maoz 
[nn](i"^ power-law best-fit). AU of the DTDs are nor- 



malized per mass formed in stars, and can be thought of 
as an absolute SN la rate as a function of Hubble time. 
We note that the observational DTDs, while currently 
the most recently-derived, have amplitudes that are at 
least a factor of a few lower than cosmic DTDs derived 
from previous studies (see Maoz et al. |153j for discus- 
sion). 

It is immediately obvious that no single progenitor 
channel can reproduce the observed rates for delay times 

< 300 Myr. The SD (red, sohd line) 'spike' at delay times 

< 200 Myr is solely due to helium-rich donors; these 
donors have relatively faster evolutionary timescales than 
their hydrogen-donor SD counterparts (hence the gap 
~ 200 - 300 Myr). The SD DTD drops off too quickly to 
follow a.t~^ power-law, but shows the expected trend at 
delay times > 400 Myr (decreasing events with increasing 
delay time). Overall, the SD rates are too low by about 
an order of magnitude to match the observed SN la rate. 

The double-detonation progenitors display a clear bi- 
modal behaviour: prompt events originate from systems 
with helium-burning star donors (green dash-dotted line; 
short evolutionary timescales) while delayed events origi- 
nate from systems where the donor is a degenerate dwarf 
(green dotted line; longer evolutionary timescales). Be- 
yond ~ 300 Myr, the rates match the observations fairly 
well, and beyond 1 Gyr the DTD follows a (steeper) 
power law ~ (see also |131j ). However, there is still 
a lot of uncertainty with respect to the explosion mech- 
anism of this channel regarding the detonation of the 
helium sheh (Sect. |IVB 2"|). 



The sub-set of DD mergers shown in Figure [2] - 
a population of violent mergers (blue dashed line) - 
are discussed in detail in Ruiter et al. jll6j (see also 
Sect. IV B 3). Once a double WD is born, it may take 



several Gyr before the stars are brought into contact. 
The dominant mechanism leading to a decrease in or- 
bital angular momentum (and hence smaller orbit) is the 
emission of gravitational waves |154j . As is expected, the 
DD mergers shown here follow the power-law shape 
fairly well, and are within the observational uncertainties 
for delay times > 300 Myr (below 300 Myr there is a dip 
in the distribution, see Sect. IV B 4 for explanation). It 



may also be possible to increase the overall rates of DD 
mergers if three-body interactions are taken into account 



IV. MODELING EXPLOSION AND 
FORMATION OF OBSERVABLES 

A. The MPA modeling pipeline 

A thorough testing of different SN la scenarios requires 
to combine several model aspects from progenitor evolu- 
tion over the explosion to the formation of the observ- 
ables. A way of addressing the progenitor problem was 
discussed in Sect. |III[ Here we present our approaches 
to modeling hydrodynamics and nucleosynthesis in the 
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FIG. 2. Coloured lines show StarTrack theoretical delay 
time distributions from 100 — 13000 Myr after star forma- 
tion assuming a 70 % binary fraction. The DTDs have been 
mass-normalized into units of SNuM on the right-hand y-axis 
for comparison with observations (SNe la {10^° Mq)"^ (100 
yr)~^). The lines have been smoothed, and small fluctuations 
are due to Monte Carlo noise. Dashed blue: DD violent WD 



mergers, which are a DD scenario sub-class (cf. Sect. IVB4I; 
solid red: SD Chandrasekhar-mass scenario including main 
sequence, giant-like and helium-burning donors; green dotted: 
double-detonation scenario (see also |131) ) with helium-rich 
WD donors; green dash-dotted: double-detonation scenario 
with helium-burning star donors. Alongside the model DTDs, 
we show the recent observationally-recovered DTD from Maoz 
et al. |153j fred squares) and the DTD from observations as 
described in Graur and Maoz [110| (straight black lines; 
best-fit). 



thermonuclear explosion and radiative transfer in the 
ejecta cloud giving rise to the observables. In the sub- 
sections which follow we first describe our combustion- 
hydrodynaniic code leafs and then the Monte-Carlo 
radiative-transfer code ARTIS. 



1. Hydrodynamic explosion models 

Type la supernova explosion models aim at following 
the hydrodynamical evolution from the ignition of ther- 
monuclear burning in a WD to homologous expansion of 
the ejecta. These models rely on the equations of hy- 
drodynamics (either the Euler equations or specific ap- 
proximations suitable for low-Mach number flows) cou- 
pled to nuclear reactions. After ignition, a combustion 
wave forms. Because of the high temperature sensitiv- 
ity of carbon fusion, it is confined to a narrow region in 
space. Seen from the scales of the WD, it can be approx- 
imated as a sharp discontinuity separating the fuel (CO 
material) from the ashes of the nuclear burning. The re- 
active Euler equations allow (in their integral form) for 
two distinct classes of discontinuous (weak) solutions that 



model the propagation of such thin combustion waves: 
subsonic deflagrations and supersonic detonations. Mi- 
croscopically, the deflagration propagation mode corre- 
sponds to a flame mediated by heat conduction, while a 
detonation is driven by a shock wave. 

The most critical aspect in modeling thermonuclear ex- 
plosions of WDs are the inherent scale problems - both 
in time and space. Although in large-scale simulations of 
supernova explosions that capture the entire WD star it 
is well justified to treat the combustion fronts as discon- 
tinuities (at least for most of the burning taking place at 
high fuel densities), this implies that the internal struc- 
ture of the fiame cannot be resolved. Thus, details of 
their mechanism and in particular the nuclear reactions 
are not represented. A time scale problem is the dis- 
crepancy between the scales of hydrodynamic flows and 
that of the nuclear reactions. In Chandrasekhar-mass 
models, the ignition of the combustion wave is preceded 
by a period of convective carbon burning (the so-called 
"simmering phase" ) that lasts for a century and is char- 
acterized by highly turbulent flows. A correct modeling 
of turbulence is also essential for deflagration phases in 
Chandrasekhar-mass explosion models. Here, Reynolds 
numbers of the order of 10^* are typically encountered 
posing another severe spatial scale problem. 

While a detailed modeling of the ignition process and 
the simmering phase in Chandrasekhar-mass explosion 
models remains challenging (but see |156L 1157] for re- 
cent efforts), there have been several attempts to over- 
come the spatial scale problem associated with combus- 
tion fronts. In large-scale supernova simulations, an 
option is to artificially broaden the combustion waves 
so that they can be represented on the computational 
grid. This flame-capturing approach was explored by 
Khokhlov 158] I159| . The main drawback is that an 
artificially-broadened flame smears out small-scale dy- 
namics. Level-set based techniques for combustion wave 
tracking are an alternative which has been introduced to 
SN la simulations by Reinecke et al. [16011161] . The un- 
derlying idea is to represent a contour (in 2D simulations) 
or a surface (in 3D simulations) by the zero-level set of 
a scalar field that is set up as a signed distance func- 
tion to the combustion wave. This scalar field is then 
evolved in an appropriate way to model the propagation 
of the burning front |160| . In this approach, the combus- 
tion wave is considered as a sharp discontinuity and no 
attempt is made to resolve its inner structure. It there- 
fore has to be augmented by a model for the propagation 
speed of the combustion wave and the energy release in 
it. For defiagrations, the laminar flame speed has been 
determined in small-scale simulations [162]. This sets the 
lower limit of the propagation speed of the effective flame 
front in large-scale supernova simulations. For most of 
the time, turbulence determines the flame speed and this 
is accounted for by a speciflc modeling approach (see be- 
low). For detonations, the propagation speed in the most 
simple case is the Chapman- Jouget speed, i.e. the sound 
speed in the ashes. This, however, does not take into 
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account the dependence on shock strength and possible 
multi-dimensional effects. In addition, at high densities 
burning to nuclear statistical equilibrium makes part of 
the process endothermic and pathological detonations are 
encountered here. For these, the microscopic mechanism 
has been studied by Sharpe jl63j and the results can be 
used to model the propagation of detonations in large- 
scale supernova simulations. 

The available computational resources limit the reac- 
tion networks employed in the hydrodynamic explosion 
simulations to only a few species. This allows to repre- 
sent the energy release to a sufficient precision to account 
for its impact on the dynamics of the explosion. The de- 
tailed nucleosynthetic yields, however, cannot be deter- 
mined directly. These are required to compare the models 
with the constraints from galactic chemical evolution and 
for setting the input models for radiative transfer simula- 
tions that rely on a spatially resolved multi-dimensional 
chemical structure of the ejecta. This problem is usually 
accounted for by a nucleosynthetic postprocessing step 
[16411165] . In the hydrodynamic explosion model, a large 
number (up to several million) of Lagrangian tracer par- 
ticles are advected with the flow. They represent fluid 
packages and for these the thermodynamic trajectories 
are recorded. These data are then used as input for a 
large reaction network that allows to reconstruct the de- 
tails of the nuclear reactions. 

Finally, the expansion of the WD and the ejecta in 
the course of the explosion pose another scale problem. 
With ejecta velocities well above 20, 000 km s^^ the mate- 
rial would quickly leave a static computational grid. Two 
approaches to overcome this problem are followed. Adap- 
tive mesh refinement allows to use large computational 
domains with finer resolution at places where physical 
processes are to be resolved. An alternative, which is im- 
plemented in our models, is to use moving computational 
meshes |166H168] . As the overall expansion of the ejecta 
is spherical to first order, a simple radial expansion of 
the computational grid provides an optimal resolution of 
the explosion physics with given computational resources 
and allows for following the hydrodynamic evolution of 
the ejecta to a relaxed state (homologous expansion) - 
a prerequisite for predicting observables with radiative 
transfer simulations. 



2. Radiative transfer models 

^From the hydrodynamic explosion simulations we ob- 
tain the velocities, densities and composition of the ex- 
plosion ejecta. These, however, are not directly com- 
parable to the observational signatures of SNe la like 
broad-band photometry, spectral time series and spec- 
tropolarimetry over a wide range of the electro-magnetic 
spectrum. For that purpose synthetic spectra and light 
curves must be obtained from radiative transfer calcula- 
tions. Since the explosion ejecta are free streaming at 
about 100s after the explosion (e.g. |166) ). the radiative 



transfer calculations can be decoupled from the hydrody- 
namic simulation assuming homologous expansion of the 
ejecta. 

The observational display of SNe la is not powered by 
the heat produced during the explosion itself. Already 
at the first observational epochs, typically a few hours to 
days after the explosion, this heat has long gone due to 
the expansion of the ejecta. Instead, the decays of ra- 
dioactive isotopes like ^®Ni and ^^Co, freshly synthesized 
during the thermonuclear burning, give rise to the emis- 
sion of a spectrum of 7-photons. These interact with the 
ejecta by Compton scattering, pair production and pho- 
toelectric absorption, thereby depositing their energy and 
reheating the ejecta [SnilSZ]- Thus, radiative transfer sim- 
ulations that aim at a direct connection between explo- 
sion models and observations have to take into account 
this energy injection and the transport of 7-photons ex- 
plicitly. A simple photospheric assumption (e.g. |169j ) is 
not enough. 

Another complication poses the peculiar chemical com- 
position of SNe la. Since their ejecta do not contain 
any hydrogen but significant amounts of iron-group ele- 
ments, the opacity in SN la is dominated by the wealth 
of lines associated with the iron-group elements (e.g. 
jl70| . figure 1), thus requiring a solution of the com- 
plicated multi-line transfer problem in expanding me- 
dia. Assuming emission from a photosphere and spherical 
symmetry, many studies have addressed this problem in 
the past either assuming pure resonance scattering (e.g. 
[US [mi [112) or pure absorption (e.g. [57]) in the fines. 
However, such an approach is too simple, since it can- 
not account for line fluorescence effects which are cru- 
cial in shaping the spectral energy distribution of SNe la 

Finally, given the complex ejecta structure of state-of- 
the art hydrodynamic explosion models a time-dependent 
3D "full-star" treatment of radiative transfer which sim- 
ulates the 7-deposition and spectrum formation in detail 
is needed. Such an approach is e.g. taken in the Monte 
Carlo radiative transfer code SEDONA of Kasen et al. [174] 
which treats line fluorescence in an approximate way sim- 
ilar to [TTDIITTH] . 

Following the methods outlined by Lucy [1751 11761 
I177j , at MPA we have developed another time-dependent 
3D Monte Carlo radiative transfer code ARTIS [1781 1179| . 
ARTIS divides the total energy available in the radioac- 
tive isotopes of a given supernova model into discrete 
energy packets. These are initially placed on a computa- 
tional grid according to the distribution of the radioactive 
isotopes and then follow the homologous expansion un- 
til they decay. Upon decay they convert to bundles of 
monochromatic 7-ray photons which propagate through 
the ejecta. ARTIS contains a detailed treatment of 7- 
ray radiative transfer [180| and accounts for interactions 
of 7-ray photons with matter by Compton scattering, 
photo-electric absorption and pair production. Assum- 
ing instantaneous thermalization of absorbed 7-ray pho- 
tons, the energy is transformed into ultraviolet-optical- 
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infrared photons enforcing statistical and thermal equi- 
librium. Using a detailed wavelength-dependent opacity 
treatment, ARTIS solves the radiative transfer problem 
self-consistently with the ionization and thermal balance 
equations. Excitation is treated approximately by as- 
suming local thermodynamic equilibrium, which is ex- 
pected to be a good approximation at least around max- 
imum light. A generalized treatment of line formation 
[T751 [T75] . including typically about 500,000 individual 
atomic line transitions |181j in the Sobolev approxima- 
tion |182| . allows for a detailed treatment of radiation- 
matter interactions including a parameter-free treatment 
of line fluorescence. Thus, depending only on the input 
model and atomic data, our radiative transfer calcula- 
tions give a maximum of predictive power for a given 
explosion model. 



B. Models for normal SNe la 

As discussed in Sect. normal SNe la can be ex- 
plained by the decay of typically 0.3 Mq to 0.9 Mgof ^'^Ni 
in the center of the ejecta which is surrounded by layers 
of intermediate mass elements, oxygen and unburnt ma- 
terial. There are several ways of constructing explosion 
models that give rise to such an ejecta structure. 



1. Chandrasekhar-mass delayed detonations 

The model of a WD exploding when approaching the 
Chandrasekhar mass is certainly the most thoroughly ex- 
plored option. Traditionally, it is associated with the 
single-degenerate progenitor model (see Sect. Ill), but 
the formation of a Chandrasekhar-mass WD due to a 
merger is not excluded. A strong argument in favor of 
the Chandrasekhar-mass model was the notion of homo- 
geneity among SNe la. This picture, however, eroded 
with the detailed observational campaigns of the past 
decade (see Sect. |ll| that clearly showed a pronounced 
diversity among these objects. By now, several sub- 
classes have been established. Therefore, it seems un- 
likely that all SNe la can be explained within a single 
progenitor/explosion model. Nonetheless, from the ex- 
plosion modeling point of view, the Chandrasekhar-mass 
scenario holds promise to explain the bulk of normal 
SNe la. There may, however, be difficulties with ex- 
plaining the rate of observed events, when this scenario 
arises exclusively in the single-degenerate channel (see 
Sect.lmi. 



One-dimensional parametrized models of explosions in 
Chandrasekhar-mass WDs have been very successful in 
reproducing normal SNe la, most notably the W7 model 
of Nomoto et al. [2]. Here, we will focus on recent de- 
velopments in simulating the nuclear burning in two or 
three spatial dimensions. There is a qualitative differ- 
ence between such models and earlier one-dimensional 
parameterizations. On the one hand, multidimensional 



approaches allow for a more realistic treatment of inher- 
ently multidimensional effects such as turbulent burning 
and asymmetries in the ignition and flame propagation. 
On the other hand, by fixing free parameters, such mod- 
els cannot easily be used to fit observations and thus 
the level of agreement with observations is usually lower. 
Thus, the assessment on the validity of the underlying 
models is more involved and interpretation is required. 

The first numerically studied explosion model, a 
prompt detonation of a Chandrasekhar-mass WD can 
be ruled out as an explanation for SNe la. Since detona- 
tions propagate at supersonic velocities with respect to 
the fuel, there is no causal contact between the energy 
release and the material ahead of the combustion wave. 
Thus, the entire star burns at the high initial densities 
(a few times 10^ g/cm'^) of a Chandrasekhar-mass WD 
in hydrostatic equilibrium. Consequently, burning pro- 
ceeds to nuclear statistical equilibrium (complete burn- 
ing) throughout most parts of the star and the ejecta 
consist almost exclusively of iron group elements (pre- 
dominantly '^'^Ni). This is in conflict with the observa- 
tional requirements for normal SNe la, see Sect. |TTj In 
order to produce less ^^Ni and a substantial amount of 
intermediate mass elements, at least parts of the burning 
must proceed at lower densities than those encountered 
in Chandrasekhar-mass WDs in hydrostatic equilibrium. 

A combustion starting out in the deflagration mode 
brings the WD out of equilibrium and pre-expands 
the fuel material. Consequently, burning partially 
takes place at lower densities than in an equilibrium 
Chandrasekhar-mass WD. This allows for the synthe- 
sis of intermediate-mass elements and reduces the ^^Ni 
yield accordingly. However, ultimately laminar deflagra- 
tion flames are too slow to catch up with the expansion 
of the star. This limits the amount of material burnt 
and thus the nuclear energy release is too low for a suc- 
cessful SN la. It has been noted early on [2] that de- 
flagrations will not propagate at their laminar speeds. 
Burning from the WD's center outward, they produce an 
inverse density stratiflcation in the gravitational field of 
the star and it is thus subject to buoyancy instabilities. 
The Ray leigh- Taylor instability and secondary shear in- 
stabilities generate strong turbulence. Driven from large 
scales, the turbulent energy cascades down to the micro- 
scopic Kolmogorov scale. Consequently, the flame inter- 
acts with turbulent eddies of various sizes. The flame 
is torn and wrinkled by these turbulent motions and this 
enlargement of the flame surface area accelerates its mean 
propagation significantly. 

Several numerical studies indicate that although tur- 
bulence is driven on large scales by buoyancy, it quickly 
becomes isotropic and follows Kolmogorov-scaling at 
smaller scales [1831 1 184] . The correct representation of 
flame-turbulence interaction is one of the key challenges 
in modeling deflagrations in WDs and thus a critical in- 
gredient in Chandrasekhar-mass models for SNe la. Sev- 
eral possibilities have been suggested to accomplish this. 
In the work discussed here, a subgrid-scale model is em- 
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ployed. It is based on a balance equation for the un- 
resolved turbulent kinetic energy. For two-dimensional 
simulations, the approach of Niemeyer and Hillebrandt 
\185[ is used while three-dimensional simulations use the 
method of Schmidt et al. |186j that does not make any 
assumptions on the scaling of turbulence. For strong tur- 
bulence, as expected for most phases of the supernova 
explosion, flame-turbulence interaction implies that on 
some sufficiently large scale (such as resolved in multi- 
dimensional simulations) the propagation speed of the ef- 
fective flame (averaged over unresolved small-scale struc- 
ture) decouples from the laminar burning speed and is set 
by the turbulent velocity fluctuations on that scale |187j . 
This is the basis for the flame model in our simulations 
that employ the level-set technique to represent the effec- 
tive deflagration front and use a subgrid-scale turbulence 
model for determining its effective propagation velocity 
(for details see also |188j ). 

The amount of burning and the energy release de- 
pend strongly on the way the flame is ignited. In 
Chandrasekhar-mass explosions a century of convec- 
tive carbon burning precedes the actual flame ignition 
(cf. Sect. IVAl). Numerical simulations of this 



phase are extremely challenging due to its long dura- 
tion and the high turbulence intensities involved (but see 
[156111891 1190 for recent attempts). At the moment, the 
geometry of flame ignition is unclear and therefore differ- 
ent possibilities are considered. If ignited in many sparks 
around the center, the WD can be unbound |1681 1191j . 
But even with a strong ignition, the asymptotic kinetic 
energy of the ejecta does not exceed ~0.6 x 10^^ erg and 
the '^^Ni production reaches at best about a third of a 
solar mass |191| . The most optimistic values for pure de- 
flagrations in Chandrasekhar-mass WDs reach the fainter 
end of normal SNe la, but they cannot account for all of 
them. Moreover, the predicted spectra show peculiarities 
that can be attributed to a chemically mixed ejecta com- 
position which is a natural consequence of the large-scale 
buoyancy instabilities in these models. Thus we conclude 
that deflagrations in Chandrasekhar-mass WDs cannot 
explain normal SNe la. They could, however, account 
for a peculiar subclass (see Sect. IV C 2). 

The only chance for Chandrasekhar-mass explosion 
models to reach the ballpark of normal SNe la is a 
detonation following the initial burning in the deflagra- 
tion mode. In contrast to a prompt detonation of a 
Chandrasekhar-mass WD in hydrostatic equilibrium, the 
pre-expansion in the deflagration phase now allows the 
detonation to burn at lower fuel densities. Although it 
still can contribute to the overall ^^Ni production, it pro- 
duces a substantial layer of intermediate mass elements 
in the outer layers of the exploding WD. One way to 
realize this is the delayed-detonation scenario |194j . in 
which a spontaneous transition of the burning front from 
deflagration to detonation occurs in a late stage of the 
explosion. This leads to a clear chemical stratification 
with iron group elements dominating the inner part of 
the ejecta while the products of a detonation in material 



of subsequently lower density lead to a stratified com- 
position in the outer layers. Here, intermediate-mass 
elements follow the iron group elements and at higher 
velocities oxygen and carbon dominate. Downdrafts of 
unburned material left behind in the turbulent and un- 
stable deflagration are now incinerated. A qualitative 
difference to earlier one-dimensional delayed-detonation 
models, however, is that stable iron group elements pro- 
duced in the high-density deflagration at the center of 
the WD do not stay there but float to larger radii due to 
buoyancy instabilities. The degree of the pre-expansion 
and thus the total ^^Ni production is determined by the 
energy release in the deflagration [SHI I195j and by the 
delay between deflagration ignition and detonation trig- 
gering. One way to vary the strength of the deflagra- 
tion is by choosing different ignition configurations (al- 
though other parameters may also affect the strength of 
the deflagration phase, see e.g. [196H199j ). Igniting vig- 
orously in many ignition sparks around the WD center 
(e.g., [16811200] ) releases more energy in the deflagration 
burning, hence achieving more pre-expansion [561 1195j . 
while a sparse and perhaps asymmetric ignition leads to a 
weak deflagration phase (e.g., [2011 1202] ). In the context 
of the delayed-detonation explosion scenario this gives 
rise to a variability of ^^Ni production which, in turn, 
leads to a range in brightnesses of the simulated events 
covering that of normal SNe la. The brightness of the 
faintest model is set by the strongest pre-expansion and 
thus by the most vigorous deflagration that is achiev- 
able. For nearly isotropic ignitions with standard WD 
setups this corresponds to a ^^Ni production in the range 
of [0.3 ... 0.4] Mq - clearly too much for subluminous 
SNe la. On the other end, weak deflagrations arising 
from asymmetric ignitions easily lead to the production 
of up to a solar mass of ^^Ni in the delayed-detonation 
scenario. Thus, in principle, this model should be able to 
reproduce the range of observed brightnesses of normal 
SNe la. 

As an example. Fig. |3] shows model NlOO 
which is ignited in 100 ignition sparks around the center. 
The ensuing deflagration (left panel) is of intermediate 
strength. The middle panel shows the deflagration front 
directly prior to the first deflagration-to-detonation tran- 
sition. The large-scale buoyancy-induced plumes of burnt 
material are clearly visible. This ~ together with shear- 
induced turbulence on smaller scales leads to the increase 
in flame surface area characteristic for the turbulent de- 
flagration. The panel on the right hand side shows a 
snapshot shortly after the flrst deflagration-to-detonation 
transition has triggered. Obviously, it is immediately fol- 
lowed by other transitions at different locations. The 
newly formed detonation waves quickly spread over the 
remaining fuel and burn out the downdrafts of fuel ma- 
terial left behind from the deflagration. Since the deto- 
nation propagates from high to low density the ash com- 
position changes from iron-group to intermediate-mass 
nuclei and, because of the supersonic propagation, there 
is no mixing, in contrast to the deflagration phase. The 
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FIG. 3. Hydrodynamic evolution of a Chandrasekhar-mass delayed detonation. Shown are a volume rendering of the density 
(orange colours) and the zero level-set of the deflagration (whitish surface) and detonation flames (blueish colours) of model 
NlOO [T921IT93] . From left to right the snapshots are taken at 0.70, 0.93 and 1.00s. 



outcome is an ejecta cloud with a stratified chemical com- 
position in the outer layers and close to 0.6 M© of ^^Ni 
at the center. The hydrodynamic evolution is followed 
with a moving-grid technique to 100 s after ignition. Af- 
ter nucleosynthetic postprocessing, the ejecta structure is 
mapped into the radiative transfer code ARTIS jl78| I179j 
to calculate synthetic observables. A sequence of spec- 
tra for this model is shown in the left panel of Fig. |4] 
Overall, the agreement between the model spectra and 
the observational reference spectra of a normal SN la 
(SN 2005cf) is reasonable. Again, we emphasize that no 
perfect match is expected in this comparison of a generic 
three-dimensional supernova model and an observation 
without any attempts of fitting. However, a more fun- 
damental shortcoming of the model is that it appears to 
be too red. This can be attributed to a flux redistri- 
bution due to stable iron group elements at rather high 
velocities - a feature that at least to a certain degree is 
characteristic for delayed-detonation models. 



A more systematic test has been presented by Kasen 
et al. |203j on the basis of a suite of two-dimensional 
models. Again, although no perfect agreement with ob- 
servational data is reached, many of the models would be 
classified as SNe la employing a tool for analyzing obser- 
vations and treating the models as actual astronomical 
data '204'. However, the brightest and most asymmetric 
explosions in the Kasen et al. I203i sample would not be 
classified as SNe la. Interestingly, in this set of models, 
the correlation between peak luminosity in the _B-band 
and the decline rate of the light curve (used to calibrate 
SNe la as distance indicators in observational cosmology, 
[To', "SS]) was found to resemble that of the observations 
p03 . Whether or not this is the case also in sets of three- 
dimensional models remains to be seen and is subject to 
forthcoming publications (Sim et al., in preparation). 



2. Sub- Chandrasekhar-mass double detonations 

The observational finding of chemically stratified 
ejecta points to a detonation propagating down the gradi- 
ent towards low densities in the outer layers of the explod- 
ing WD. As discussed above, for Chandrasekhar-mass 
WDs this is only compatible with a configuration that 
is out of hydrostatic equilibrium. An alternative to this 
mechanism is a detonation in a sub-Chandrasekhar mass 
WD. Pure detonations in CO WDs with masses between 
0.81 Mq and 1.15 M© have been tested by Sim et al. [207] 
(see also Shigeyama et al. [208] ) and yield ^^Ni masses in 
the range of [0.01 . . . 0.81] M©. According to the model 
sequence of Sim et al. |207j , a standard normal SN la with 
^ 0.6 Mq of '^^Ni is expected to result from a detonation 
in a WD of about 1.1 M©. The observables predicted from 
these models roughly match the data from normal SNe la 
and their B-band light curves seem to follow the width- 
luminosity relation [207! • Thus, detonations in WDs with 
masses well below the Chandrasekhar-limit hold promise 
for explaining normal SNe la. The question is how a 
detonation in such an object can be triggered. Here we 
discuss one possibility arising from a detonation in an 
accreted He shell on top of the WD. Another possibility 
- due to the merger of two WDs - will be presented in 
the next section. 

The idea of double detonations in sub-Chandrasekhar 
mass WDs has been discussed extensively in the 1990s by 
Woosley and Weaver [209] , Livne and Arnett [210J , Benz 
[m] . Livne [US, Garci'a-Senz et al. [HS]- A CO WD ac- 
cretes helium from a companion star (either a helium star 
or a helium WD). When the accreted He layer becomes 
sufhciently massive, compressional heating is thought to 
lead to a detonation in the He material (see, however, 
|214j for an alternative mechanism based on instabili- 
ties in the accretion process). This detonation sweeps 
around the CO core and burns the He to heavier ele- 
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FIG. 4. Synthetic spectra of different kinds of explosion models for normal SNe la. From left to right the panels show (i) the 
delayed-detonation model NlOO [192| . (ii) model 3m of the sub-Chandrasekhar-mass double detonations presented by Kromer 
et al. |2U5j . and (iii) a double-degenerate merger of two WDs with 1.1 and 0.9 Mq [2U6) . For comparison, we show observed 
spectra of the "golden-standard" normal SN la 2005cf for corresponding epochs [40] (data in red). 



mcnts. At the same time a shock wave propagates into 
the core. This shock may trigger a secondary detonation 
close to the interface between the CO core and the He 
shell ("edge- lit detonation"), or when reaching the cen- 
ter of the core. The secondary detonation incinerates the 
entire WD and leads to its successful disruption in a ther- 
monuclear supernova. The question, however, is whether 
the event would really look like a SN la. Although for 
sufficiently massive CO cores enough ^®Ni can be pro- 
duced to power a normal SN la, problems arise from the 
burning products of the He shell. In the models of the 
1990s, a rather massive He shell ~ about [0.1 . . . 0.2] Mq - 
was thought to be necessary to trigger a detonation and 
to drive a sufficiently strong shock wave for initiating 
a secondary detonation in the core. In such massive He 
shells, a detonation produces a significant fraction of iron 
group elements (including additional '^^Ni). These affect 
the radiative transfer and the predicted observables are 
at odds with the actual observations |215H218] . 

Recently, however, Bildsten et al. |220j and Shen and 
Bildsten 221^ pointed out that in AM CVn systems 
rather low masses of accreted He on top of a CO WD can 
develop dynamical burning, possibly in the detonation 
mode. The work by Fink et al. |219L I222j demonstrated 
that a core ignition is very robust due to spherical shock 
convergence near the center of the WD which leads to a 
geometrical shock amplification. Neither asymmetric ig- 
nition geometries [222] nor low He shell masses prevent a 
secondary core detonation once the He shell successfully 
triggers a detonation [219] . The low He shell mass sig- 
nificantly reduces the observationally disfavored effects 



of iron group elements in the outer layers of the ejecta 
[205] . In addition to lowering the total mass of the He 
shell, it also reduces the density at which He detonates 
thus leading to predominantly incomplete burning. Con- 
sequently, the outer layers of the ejecta in the models of 
Fink et al. 219.^ contain virtually no ^^Ni and only low 
amounts of other iron group elements. Another effect 
that distinguishes the models of Fink et al. |219[ from 
many earlier calculations and also from the recent mod- 
els of Woosley and Kasen |218j is the multi-dimensional 
treatment of the He shell detonation. Sweeping around 
the CO core, it propagates laterally and allows for sig- 
nificant post-shock expansion - an effect that is not cap- 
tured in spherically symmetric models. This adds to the 
less complete burning observed in the simulations of Fink 
et al. |219j (see also 223 ). According to Kromer et al. 
[205] the reduced yields of heavy elements from the He 
detonation have significant impact on the predicted ob- 
servables. Although the colors are too red to match the 
observations perfectly, the range of brightnesses and rise 
and decline rates of normal SNe la is covered by the mod- 
els. 

Kromer et al. [205] also point out that the details of 
spectra and the colors are very sensitive to the thermal 
and chemical conditions in the detonating He shell. In 
particular, they find that the degree of burning in the 
shell material (and thus its final composition) can be af- 
fected by the initial composition of the He shell. Since 
the time-scale for a-captures behind the detonation shock 
front is significantly shorter than that of triple-a reac- 
tions, a ^^C admixture in the He shell due to previous 
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FIG. 5. Hydrodynamic evolution of a sub-Chandrasekhar- 
mass double detonation for a helium shell of 0.084 M© on top 
of a 0.920 M0 CO WD (from top left to bottom right; model 
2 of Fink et al. [2T9| '). The density structure of the WD is 
colour coded in red {pi in units of lO^gcm"^). The dashed 
blue lines indicate the border of the helium shell. The solid 
blue and magenta lines show the helium and CO detonation 
flames, respectively. 



hydrostatic burning or dredge-up of core material [221] 
can limit the a-chain before reaching nuclear statistical 
equilibrium. In an exploratory model Kromer et al. j205j 
homogeneously polluted a He shell with 34% (by mass) of 
^^C (their model 3m) and showed that such a model pro- 
duces light curves and spectra that are in good agreement 
with those of normal SNe la. In particular, this model 
is no longer too red at maximum light. However, one 
caveat remains: the model still produces a non-negligible 
amount of Ti in the outer layers leading to the formation 
of a Ti II absorption trough between 4000 A and 4400 A 
which is not observed in normal SNe la (see also middle 
panel of Fi g. [4|) b ut in subluminous 91bg-like supernovae 



only (Sect. II B I. Whether these differences can be re- 
solved as well, remains to be seen in future studies that 
more fully explore the influence of the initial composi- 
tion of the helium shell and different ignition geometries. 
Also the strong sensitivity of the radiative transfer to 
tiny amounts of particular elements requires a better de- 
scription of nuclear reaction rates and continued study 
of the radiative transfer processes (and atomic data) in 
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FIG. 6. Hydrodynamic evolution and subsequent thermonu- 
clear explosion of a merger of a pair of 1.1 and 0.9 M© WDs 
[206]. Initially the WDs orbit each other with a period of 
~ 35 s. After a few orbits the secondary is tidally disrupted 
and collides with the more massive primary reaching densi- 
ties and temperatures sufficient to ignite a detonation at 610 s 
(black cross). At 612s the detonation front (black line) has 
burned almost the complete object. Color-coded is the log- 
arithm of the density. Note that the last two panels have a 
different color scale ranging from 10"* g cm~^ to 10^ g cm~^ 
'g cm~^, respectively. 



order to quantify more fully the systematic uncertainties 
which arise due to the complexity of spectrum formation 
in supernovae. 



3. Violent mergers 

Another external trigger to ignite a sub- 
Chandrasekhar-mass WD is the violent merger of 
two CO WDs. Although the total mass of the merging 
system usually exceeds the Chandrasckhar mass, both 
components are below this mass limit. The so-called 
violent merger model [224] starts with massive WDs 
[M > 0.9 M0 ) with a mass ratio close to unity. This 
scenario results from a subset of double-degenerate 
progenitor models. Other configurations of merging 
WDs may avoid thermonuclear explosions (e.g., |225j ) 
and instead lead to the formation of a neutron star by 
gravitational collapse |144| . 

For mass ratios close to unity, however, the merg- 
ers proceed dynamically and can be followed in hydro- 
dynamic simulations. Pakmor et al. |206j studied the 
merger of a 1.1 Mq "primary" WD with a 0.9 Mq "sec- 
ondary" WD. The inspiral and merger, as followed with 
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the SPH code Gadget [226] in its modification for stel- 
lar astrophysical problems |227j is shown in Fig. [6j In the 
last few orbits before the actual merger, tidal interaction 
strongly deforms the secondary and it finally plunges into 
the primary WD (snapshots for t > 600 s in Fig. [6]). This 
violent merger leads to the formation of a hot spot where 
the two masses collide (marked by a black + in the snap- 
shot for t — 610 s). Here, thermodynamic conditions are 
suitable for triggering a detonation. 



After mapping into our grid-based supernova explo- 
sion code, the detonation (indicated by a black contour 
in Fi g. [6| is followed with the level-set technique (see 
Sect. |IV A 1|). It incinerates the merged object almost 
An important point to notice is that the 



completely. 

primary WD is nearly unaffected by the merger. There- 
fore the burning takes place at the low densities typi- 
cally found in sub-Chandrasekhar-mass WDs. Only the 
primary possesses material at sufficiently high densities 
to synthesize iron group elements while the secondary 
mostly burns to oxygen. With the mass of the primary 
chosen to 1.1 Mq, a moderate ^®Ni mass production is ex- 
pected according to |207j and indeed 0.64 of ^^Ni are 
found in the presented simulation. Significant amounts of 
the total 2 of material involved in this merger burn to 
intermediate-mass elements (0.6 M©) and the ejecta con- 
tain 0.47 M0 of oxygen. Only 0.09 M0 of carbon remain 
in the ejecta. 



Thus, despite the large total mass of the exploding ob- 
ject, the angle-averaged and line-of-sight dependent light 
curves of this merger compare very favorably to that of 
normal SNe la. With a peak brightness of —19.6, —19.0, 
and —19.2 in the U, B, and V bands, respectively, and 
a B-band light curve dechne rate of Ami5(B) = 0.95, 
the model predictions are well in the range of those ob- 
served for normal SNe la [57] • Moreover, the spectral 
evolution of the model (see Fig.|4]) reproduces the overall 
spectral shape and the velocity-shifts of most of the line 
features remarkably well and shows most of the charac- 
teristic features of SNe la, particularly the defining Sill 
doublet at AA6347,6371 but also the weaker Sin features 
at AA5958,5979 and AA4128,4131. Other prominent fea- 
tures are the Call H and K absorptions, the Mgii triplet 
at A4481, the Sii W-feature at ~5400A and, in the red 
tail of the spectrum, the O i triplet AA7772, 7774,7775 and 
the Can NIR triplet. However, with a B-band rise time 
of 20.8 d the pre-maximum light curve evolution of this 
merger model is relatively long compared to that of nor- 
mal SNe la. Hayden et al. [Hj, for example, find an 
average B-band rise time of 17.4 d in the SDSS-II SN 
sample (but see also Conley et al. |33i who find a value of 
19.58 d for a low red-shift sample from the SNLS). This 
could indicate that the total ejecta mass in this model 
is somewhat too large. Future studies exploring the pa- 
rameter space of violent WD mergers in more detail will 
show if this explanation is right. 



4- Critical assessment 

Although the models presented above produce success- 
ful explosions that overall compare favorably to the ob- 
servations of normal SNe la, there remain uncertainties 
in the modeling of the explosion physics. In all mod- 
els this refers to the initiation of the burning and the 
formation of detonations. This is not too surprising as 
these processes work on scales that cannot be resolved 
in our multi-dimensional supernova simulations. More- 
over, the ignition of deflagrations and detonations are 
complex physical phenomena and their microphysics not 
completely understood, even for terrestrial combustion. 
Despite our attempt to model the explosion physics as 
parameter-free as possible, we are thus left with the fol- 
lowing critical points in the three classes of models: 

• Delayed detonations in Chandrasekhar-mass 
WDs hinge on the possibility of deflagration-to- 
detonation transitions to occur in WD combustion. 
Although some recent studies (e.g., |228H231] ) 
indicate that this may indeed be the case, it is 
difficult to definitely decide on its realization in 
SNe la. The other major uncertainty in this model 
is the way the defiagration ignites. A strongly 
asymmetric ignition leads to extremely bright 
events in the context of delayed detonations. 
Reaching the low-luminosity end of the normal 
SNe la requires to limit the ^^Ni production to the 
yield of the strongest pure deflagrations. These 
result from quasi-isotropic or central ignitions of 
the WD - a scenario that is currently not favored 
by ignition simulations [156 but may arise due to 
slight rotation in the ignition phase [190] . 

• Double detonations in sub-Chandrasekhar mass 
WDs require the initiations of two detonations. 
While the secondary detonation in the CO core 
seems to be virtually unavoidable |219l I222| . the 
initial detonation in the He shell is not established 
beyond doubt - in particular for the case of low 
He-shell masses. 

• Violent mergers of two WDs rely on the trigger- 
ing of a detonation at the encounter of the two 
stars. Although the simulations of Pakmor et al. 
[2061 12241 1232] indicate that this is possible, the 
mechanism still awaits a firm proof. 

A better understanding of the microphysics of ther- 
monuclear combustion is thus required to overcome these 
uncertainties and to assess the models purely from the 
plausibility of their explosion mechanism. This is a chal- 
lenging task and a convincing result is not expected in 
the short term. There are, however, alternative ways to 
judge the potential of different explosion scenarios to ac- 
count for the majority of normal SNe la. 

One possibility, which we followed here, is to per- 
form supernova simulations under the assumption that 
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the uncertain mechanisms in the modehng proceed in 
a favorable way and to compare the outcome with ob- 
servations. Our results indicate that all of the models 
considered here are able to cover the range of explo- 
sion energies and brightnesses of normal SNe la and to 
first order reproduce their light curves and spectra rel- 
atively well. The comparisons of our synthetic spectral 
time series and SN 2005cf in Fig. |4] demonstrate this suc- 
cess, but also show that in detail there are shortcom- 
ings in each of the models as discussed in the previous 
sections. In the Chandrasekhar-mass delayed-detonation 
and sub-Chandrasekhar-mass double-detonation models 
the blue-shift of the characteristic Si ii is too large com- 
pared to SN 2005cf, indicating slightly too high ejecta 
velocities. This potentially can be cured by more real- 
istic progenitor models with carbon-depleted cores |192j . 
Moreover, the models are too red compared to the ob- 
servations. This is most pronounced for the delayed- 
detonation Chandrasekhar-mass model but also found in 
the double-detonation sub-Chandrasekhar-mass and the 
violent merger model to some degree. 

Although the involved masses and the explosion 
physics of the models shown in Fig. |4] differ significantly, 
at the current precision of the models it is difficult to dis- 
tinguish them by means of maximum-light optical spec- 
tra only. This degeneracy prevents favoring one model 
over others. This could imply that all channels con- 
tribute to normal SNe la (possibly with different real- 
ization frequencies) or other ways of discriminating them 
have to be found. Promising for this task seem observa- 
tions in the ultraviolet (e.g., |105l 1233] ) and in the near- 
infrared bands (e.g., [234]) but also late time observations 
(e.g., (192i i235j ). spectropolarimetry (e.g. . 12361 1237j ) or 
gamma-ray observables (e.g., |180[ 12381 1239| ). For these, 
either theoretical models have yet to be developed for 
modern multi-dimensional SN la simulations, or data has 
to be acquired. 

Other possibilities to discriminate different explosion 
scenarios are the search for signatures of the progenitor 
system in nearby SNe la (see, e.g., [HI [Ml L52 , 53, 240, 
I241| for observational constraints and, e.g., [242n245j for 
theoretical predictions) or in supernova remnants (e.g.. 

Finally, any model scenario that is claimed to account 
for a large fraction of SNe la must be able to explain ob- 
servational trends like e.g. the observed delay-time and 
brightness distribution of SNe la. By combining the syn- 
thetic observables from our explosion models with studies 
of the realization frequency of the supposed progenitor 
systems, as discussed in Sect. we can thus put addi- 
tional constraints on the different explosion scenarios. 

Recently, we have used this approach to investigate the 
prospect of the violent merger scenario in more detail. As 
was noted in Sect. |IV B 3[ the secondary WD in a violent 
merger - while consumed in the explosion - does not 
contribute to synthesizing ^^Ni. Thus, it is the primary 
WD (more specifically its mass) that simply determines 
the peak luminosity of a SN la in the violent merger 



model. 

For a realistic estimate of primary WD masses in 
would-be merging WD pairs, we took the distribution 
of primary WD masses of all merging CO WDs from 
the binary evolution population synthesis calculations of 
Ruiter et al. |131j (their standard model). Using this 
mass distribution, a relationship between the (primary, 
sub-Chandrasekhar mass) WD and its corresponding SN 
bolometric peak brightness (mwD ^ -'^^boi) was derived 
using the technique as described in Sim et al. {2071 . 

One critical question is the realization of a WD merger 
itself: e.g. what is the critical mass ratio for which mass 
transfer will be dynamically unstable (and lead to a 
merger) when the larger WD fills its Roche-lobe? The an- 
swer to this question is not straightforward, and much un- 
certainty exists in the modeling of mass transfer in close 
binaries |250| . A trend that was found in the merger sim- 
ulations of Pakmor et al. [206] (see also |25lJ) is that the 
critical mass ratio qc likely decreases with larger primary 
masses. We constructed a relationship f |116j . equation 1) 
that follows this trend to evaluate whether a given double 
WD could produce a merger that is sufficiently violent. 
Additionally, for the violent WD mergers we limited the 
primary mass to be above 0.8 Mq since primaries less 
massive than this are expected to barely produce even 
0.01 Mq of s^Ni ([203, table 1). 

Fig. [7] shows four model (bolometric) peak brightness 
distributions for a range of ^c-cuts. In grey scale we over 
plot the (scaled-up) observational luminosity distribution 
of SNe la from Li et al. |33j. Regardless of the assumed 
Qc-cut, our theoretical brightness distributions do a fairly 
good job in covering the range and matching the shape 
of the observed SN la brightness distribution. Such good 
agreement indicates that merging WDs which explode via 
the violent merger mechanism could be dominant SN la 
progenitors, driving the shape of the underlying bright- 
ness distribution. 

In Fig. [8] we show the primary mass in the same model 
for violent WD mergers as a function of delay time. The 
darkest hexagons represent the regions of highest density 
for a given cell. It is clear from the plot that mergers 
hosting the most massive primaries (> 1.3 M©) tend to 
merge at prompt (< 500 Myr) delay times. (We note 
that in single star evolution, CO WDs would not achieve 
such high masses and a WD of mass ^1.3 Mq would 
be composed mostly of oxygen and neon. However, in 
binary evolution, such masses are allowed for CO WDs, 
in particular if the CO WD accretes mass after it has 
formed (see e.g. 116 , fig. 2)). Since in the violent merger 
model the SN la luminosity is determined by the primary 
WD mass, this means that we would expect the bright- 
est SNe to be found amongst very young stellar pop- 
ulations; a trend which is confounded by observations 
[252] . One aspect of the binary evolution model which 
remains to be confirmed, is whether or not primary WDs 
are able to efficiently accrete on the order of 0.2 Mq 
from a slightly-evolved helium star companion. Such a 
mass transfer phase was found to be critical in producing 
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FIG. 7. Brightness distribution of violent WD mergers [see 
11161 for details]. Black solid histogram shows all CO WD 
mergers from population synthesis, while coloured histogram 
lines show the brightness distributions when more stringent 
mass ratio constraints are assumed. Grey scale shows the 
observational peak brightness distribution of 74 SNe la from 
the volume-limited sample of Li et al. [33]; observations are 
scaled up to enable comparison with the distribution shapes 
from our models. 



a large number of primary WD masses that yield peak 
explosion brightnesses around —19 mag (see |116j ). Ad- 
ditionally, we note that there exists another population 
of violent mergers with very short (< 100 Myr; ultra- 
prompt) delay times. These systems undergo two com- 
mon envelope phases, whereby the secondary star loses 
its envelope twice. 



Peculiar SNe 



As discussed in Sect. |IIB[ as of today several pecu- 
liar sub-classes of SNe la have been found in addition to 
the bulk of spectroscopically normal SNe la which follow 
the Phillips relation. Here, we discuss possible explosion 
models for a few of those peculiar sub-classes. 



1. 1991bg-Uke SNe 

1991bg-like SNe are subluminous with respect to the 
Phillips relation and peak at about -17 mag, indicating 
that only a rather low ^^Ni mass of about 0.1 M© was syn- 
thesized during the explosion. Moreover, a spectral anal- 
ysis of SN 2005bl, a well-observed proto-typical 1991bg- 
like SN, has shown that both iron group elements and 
silicon are present over a wide range of radii extending 
down to very low expansion velocities. This indicates the 
presence of incomplete Si burning over a wide velocity 
range in these explosions as it may occur in detonations 
at low densities. 
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FIG. 8. Delay time distribution for primary WD masses at 
time of merger for one of our violent merger model popula- 
tions (cf. blue histogram in Fig. [7|. We show delay times only 
from — 2000 Myr so that the characteristic primary masses 
at prompt delay times are clearly visible. The most massive 
primaries tend to merge at delay times < 400 Myr, however 
there also exists a distinct population of 'ultra-prompt' merg- 
ers with less-massive primaries with delay times < 100 Myr 
(see text). 



In the violent merger scenario (see also Sect. IV B 3 



]) such a burning is possible for a primary WD with a 
sufficiently shallow density profile. Following the inspiral 
of a pair of 0.89 Mq WDs with the SPH code GADGET 
and using our full modeling pipeline Pakmor et al. |224| 
have shown that such a configuration produces about the 
right amount of ^^Ni although a total mass of 1.8 Mq 
is involved in the merger. Moreover, their simulation 
can reproduce the observed spectra and light curves of 
1991bg-like SNe (see Figure [9]) and accounts for most of 
their peculiar features. 

Pakmor et al. [232) find that WD binaries with a 
primary mass of Mi ~ 0.9 Mq and mass ratios q — 
M2/M1 > 0.8 evolve similarly, thus confirming a robust 
ignition of ^0.9 Mq violent WD mergers which makes 
them promising candidates for 1991bg-like SNe given 
the good agreement of synthetic observables in our pilot 
study. Arguing that primary WDs with a lower mass will 
not detonate due to their lower densities and using popu- 
lation synthesis calculations of Ruiter et al. [126 , Pakmor 
et al. |224| also estimated the rate of binary mergers that 
met the necessary criteria to satisfy their model. It was 
found that such mergers may contribute on the order of 
2 — 11 % to the total SN la rate, which is not too far off 
from the observationally derived rate for 1991bg-likes of 
15% [33]. Moreover, if the WD binaries undergo only one 
common envelope phase and/or begin their evolution on 
the ZAMS with wide orbital separations their model also 
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FIG. 9. Angle-averaged synthetic light curves of a merger of two 0.89 M0 WDs (black). To indicate the spread due to different 
viewing angles, the gray lines show light curves along four different lines-of-sight. These have been selected from 100 equally 
sized solid-angle bins such that they represent the full range of the scatter. For comparison observed photometry for normal 
(blueish colours, Krisciunas et al. [253], Pastorello et al. [25411255] ) and sub-luminous 1991bg-like SNe (red, Taubenberger et al. 
[59] and references therein) is shown. 



prefers old (> 1 Gyr) stellar populations as indicated by 
observations of 1991bg-like SNe. 



2. 2002cx-like SNe 

One of the most peculiar sub-classes of SNe la are 
explosions similar to SN 2002cx [75]. Those events are 
under-luminous with respect to the Phillips relation and 
their NIR light curves do not show secondary maxima. 
Moreover, their spectra are characterized by very low ex- 
pansion velocities compared to normal SNe la and show 
signs of strongly mixed ejecta. While explosion mod- 
els involving a detonation are not able to explain such 
an ejecta structure (e.g., [1931 1207j ). turbulent deflagra- 
tions in Chandrasekhar-mass WDs naturally predict such 
a strong mixing and low kinetic energies [256, 257J. 

Given our ignorance of the exact ignition configu- 
ration of Chandrasekhar-mass WDs (see Sect. IVB 1| , 
we have recently performed a systematic study of 3D 
full-star explosion simulations of pure deflagrations in 
Chandrasekhar-mass WDs (Fink et al., in preparation) 
for different ignition setups. Depending on the strength 
of the ignition which is parametrized by a varying number 
of ignition sparks to seed unstable burning modes we ob- 



tain ^^Ni masses between 0.035 and 0.38 M©. Moreover, 
we find that only strong ignitions release enough energy 
during the burning to unbind the progenitor WD com- 
pletely. Asymmetric, weak ignition setups, in contrast, 
lead to a one sided deflagration plume which fragments 
due to Rayleigh- Taylor and Kelvin-Helmholtz instabil- 
ities and finally wraps around the still unburned WD 
core when it comes close to the surface (Figure 10 see 
also I 201[ I258[ ). However, even deflagrations which fail 
to unbind the complete WD accelerate parts of their ex- 
plosion ashes to escape velocity and eject this material 
into their surroundings. 



Using a million Lagrangian tracer particles, we deter- 
mined the detailed chemical composition of our simula- 
tions from a post-processing calculation with our 384- 
isotopes nuclear network [1(31 1165| and mapped the re- 
sulting ejecta structure into our radiative transfer code 
ARTIS |178[ I179[ . While the obtained synthetic observ- 
ables for strong deflagrations, which completely unbind 
the progenitor WD, do not match the display of observed 
SNe, deflagrations which leave behind a bound remnant 
closely resemble the observed properties of 2002cx like 
SNe (Figure 11 for details see Kromer et al. [259| and 
Jordan et al 
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FIG. 10. Hydrodynamic evolution of an asymmetrically ignited deflagration in a Chandrasekhar-mass WD. Particularly we 
show model N5def of Fink et al. (in preparation) which ejects about 0.37 M© of which ~ 0.16 M0 are ^^Ni. Shown is a volume 
rendering of the mean atomic number (colour bar) from which we carved out a wedge to allow a view into the core. At 0.75 s 
after the explosion a one-sided deflagration plume rises towards the WD surface which fragments due to buoyancy instabilities, 
(ii) At 1.5 s the expansion of the WD quenches the burning and the explosion ashes wrap around the unburned core, (iii) 
Finally, at 100 s the unburned core is completely engulfed by the explosion ashes which are accelerated to escape velocity. 




4000 6000 8000 

Rest wavelength In A 

FIG. 11. Synthetic spectra of the asymmetrically ignited de- 
flagration model N5def which leaves behind a bound remnant. 
The spectral evolution is remarkably similar to SN 2Q05hk 
|60| , a proto-typical 2002cx-like SN. For comparison we show 
also spectra of SN 2005cf |40j as an example for a normal 
SN la. 



3. Superluminous or "super- Chandra" SNe 

Recently observations revealed a new class of superlu- 
minous SNe la (e.g.,[STJ[5n]) with total ejecta masses sig- 
nificantly larger than the canonical Chandrasekhar mass 
of 1.4 Mq. SN 2009dc, which is one of those objects, re- 
quires even a ^®Ni mass larger than the Chandrasekhar 
mass if its peak luminosity was solely powered by ra- 



dioactive decay of ^^Ni and its daughter nuclei [6T| [95]. 
One model which was proposed for these peculiar objects 
by Howell et al. [SO] is that of exploding rapidly rotating 
WDs which stay stable well above 1.4 Modue to centrifu- 
gal forces |261j . 

Several authors have studied prompt detonations |2621 
I263j and turbulent deflagrations |264j in such differen- 
tially rotating WDs. Here we report on a delayed deto- 
nation in a differentially rotating WD of 2 (see Fig- 
Fink et al 



ure 
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in preparation). Compared to a de- 
layed detonation in a non-rotating WD, the initial defla- 
gration propagates preferentially along the rotation axis 
since angular momentum conservation and weaker gradi- 
ents in the effective potential inhibit the growth of flame 
instabilities in lateral directions. A similar effect was al- 
ready found by Pfannes et al. [264] for pure deflagrations 



in such an object. As a consequence not much energy 
is released during the deflagration phase. Therefore, the 
WD does not expand strongly before the deflagration- 
to-detonation transition leaving a large amount of fuel 
at high densities which the ensuing detonation efficiently 
burns to nuclear statistic equilibrium. 

Yielding a total ^^Ni mass of 1.45 our simulation 
gives rise to a bright explosion which in principle qualifies 
the model as an explanation of super-luminous SNe la. 
However, the observationally derived ejecta structure of 
those objects [SS] does not match our explosion. This 
is also reflected by the synthetic observables from our 
model which do not match the observed spectra of super- 
luminous SNe la (see Figure 13 ). In particular absorption 



features of intermediate-mass elements, such as Si and S, 
are significantly blue-shifted with respect to the observed 
spectra, thus indicating that these elements are located 
at too large velocities in our model. Moreover, we cannot 
reproduce the characteristic C features of super-luminous 
SNe la in our model since the detonation burns almost 
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FIG. 12. Delayed detonation in a differentially rotating WD 
of 2M0. Sliown is a snapshot of the flame evolution at 0.9 s 
after the explosion. It is clearly visible that the propaga- 
tion of the wrinkled deflagration flame (reddish surface) is 
inhibited in lateral directions by the rotation and propagates 
predominantly along the rotation axis. At several location 
a deflagration-to-detonation transition occurred and detona- 
tion flames started to spread (whitish surface). The donut 
structure of the differentially rotating WD is indicated by the 
blueish volume rendering of its density. 
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FIG. 13. Synthetic maximum light spectrum of a delayed 
detonation in a differentially rotating WD of 2M0. The char- 
acteristic absorption features of Si and S are blue-shifted with 
respect to the observed spectrum of SN 2009dc [61] , a proto- 
typical super- luminous SN la. The binding energy of the WD 
is not large enough to compensate for the huge energy release 
due to nuclear burning leading to a too high kinetic energy. 



all the fuel. 

As an alternative explanation for super-luminous 
SNe la Howell et al. proposed the merger of two 
massive CO WDs. However, at least in the violent WD 
merger scenario |224j . this seems to be unlikely. In this 
model the produced ^^Ni mass depends only on the mass 
of the primary WD. Since exploding CO WDs in the vio- 
lent merger model usually have masses well below 1.3 



[116] , this essentially limits the achievable ^^Ni mass in a 
violent merger to ~ 1 . For an alternative explanation 
of super-luminous SNe la in an interaction scenario see 
Hachinger et al. W]. 



V. SUMMARY AND CONCLUSIONS 

In this article we have reviewed some of the recent 
work on SNe la done by the MPA-Garching group. Most 
of this work was motivated by the fact that this class of 
stellar explosions is not as homogeneous as it appeared 
to be in the past. In fact, new detailed observations of 
many nearby events as well as results from recent super- 
nova surveys seem to indicate that there is not a single 
progenitor channel but that several distinctively different 
channels are more likely. 

Therefore we have started a new effort to simulate 
not only single-degenerate Chandrasekhar-mass explo- 
sions but also sub-Chandrasekhar mass models and (vi- 
olent) double-degenerate mergers, to compute synthetic 
light curves and spectra from these models, and to com- 
pare their predictions with data. Moreover, we have 
made an attempt to compute rates and delay times of 
the different progenitor classes from binary-population 
synthesis models. The main results of this program were 
presented in the previous sections. 

We have demonstrated that some of the models are 
able to reproduce light curves and spectra of 'normal' 
SNe la amazingly well, given the fact, that these mod- 
els have almost no tunable (non-physical) parameters. 
The agreement is not perfect but, given the uncertainties 
still present in the models (initial conditions, combustion 
physics, radiative transfer, ...) this is not so surprising. 
Also, it may be better not to compare an individual su- 
pernova with a particular realization of a special group 
of models but try to reproduce 'generic' features of a 
full class of objects instead. This will become possible 
in the future once extended grids of models have been 
computed. 

The bad news is that rather different explosion mod- 
els reproduce the data equally well (or not so well). This 
'degeneracy' can be understood from the fact that mainly 
the abundances and distribution of radioactive ^^Ni and 
intermediate-mass nuclei determine the observed prop- 
erties of thermonuclear supernovae, and they are not 
too different for delayed-detonation Chandrasekhar-mass 
models, for sub-Chandrasekhar-mass explosions, or for 
violent mergers. These classes of models differ mainly in 
their ejecta masses and the amount of unburnt carbon 
and oxygen. But since the opacity of C and O is low this 
has, in general, little effect on the light curves and spec- 
tra. In the future, strong arguments in favor of one or 
the other progenitor channel may come from constraints 
on the rates and delay times. Our population synthesis 
models support the double-degenerate scenario, but this 
is still controversial. Other constraints come from direct 
observations, such as the presence or absence of circum- 
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stellar gas, the non-detection of the progenitor star or its 
companion, and so on. As it stands, the results of such 
studies are conflicting, but could best be explained by 
more than one progenitor channel. 

As far as some of the peculiar SNe la are concerned 
models appear to be more conclusive. We have shown 
that SN 2002cx-like supernovae can be explained well 
by pure deflagrations of Chandrasekhar-mass WDs that 
leave behind a bound WD. SN 1991bg-likc events, in con- 
trast, can be explained by a violent merger of two WDs 
of almost equal mass around 0.9 Mq. Finally, in our sim- 
ulations we did not find an explanation of the superlumi- 
nous SNe la. Neither the merger of two massive WDs nor 
the explosion of a rapidly-rotating super-Chandrasekhar- 
mass WD can reproduce the high luminosity and rather 
normal expansion velocity at the same time. Here, one 
might speculate that not all the luminosity comes from 
radioactive decay. 
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